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Abstract 

We review some recent results that express or rely on the locality properties of the dynamics 
of quantum spin systems. In particular, we present a slightly sharper version of the recently 
obtained Lieb-Robinson bound on the group velocity for such systems on a large class of metric 
graphs. Using this bound we provide expressions of the quasi-locality of the dynamics in various 
forms, present a proof of the Exponential Clustering Theorem, and discuss a multi-dimensional 
Lieb-Schultz-Mattis Theorem. 



1 Introduction 

Locality is a fundamental property of all current physical theories. Sets of observables can be 
associated with points or bounded regions in space or space-time and a relativistic dynamics will 
preserve this structure [8j. How the locality property manifests itself mathematically in important 
situations continues to be an active topic of investigation |13j . 

There is a wide range of important physical systems, however, which we prefer to describe by 
very effective non-relativistic quantum theories with Hamiltonian dynamics. Even if the Hamil- 
tonian has only finite-range interactions, the dynamics it generates generally does not preserve 
locality, i.e., there is no strict equivalent to the finite speed of light. However, locality still holds 
in an approximate sense, and there is an associated finite velocity, which is sometimes referred to 
as the group velocity. We call it the Lieb-Robinson velocity since Lieb and Robinson were the 
first to prove its existence and to obtain a bound for it [16j. They proved that to a high degree 
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of accuracy locality is preserved by quantum spin dynamics in the sense that any local observable 
evolved for a time t > remains localized in a region of space with diameter proportional to t, up 
to an arbitrarily small correction. This also means that spatial correlations between observables 
separated by a distance d cannot be established faster than a time of order d. 

The fundamental issue of locality may be sufficient motivation to extend the Lieb-Robinson 
bounds to more general situations, but there are other good reasons to try to generalize their result 
and to improve the estimates they obtained. As we will discuss below, locality, or the approximate 
locality of the dynamics, has been shown to be responsible for a considerable number of other 
important properties relevant for models of many-body systems. In many situations, however, the 
implications of locality have yet to be fully explored. 

We will begin this note by presenting a short proof of the new Lieb-Robinson bounds obtained 
successively in [21], [12], and [19]. This improved result we give below sharpens the bounds previ- 
ously obtained in that the prefactor now only grows as the smallest surface area of the supports of 
the local observables. We do this in Section [2j An application where this surface area dependence, 
rather than volume dependence, is important can be found in [6]. 

In Section [3l we present two perspectives on how Lieb-Robinson bounds may be used to provide 
explicit estimates on the local structure of the time evolution. As a consequence, one easily derives 
bounds on, for example, multiple commutators and the rate at which spatial correlations can be 
established in normalized product states. 

Section [5] discusses the so-called Exponential Clustering Theorem. In the relativistic context 
it has been known for a long time that a gap in the spectrum above the vacuum state implies 
exponential decay of spatial correlations in that state [2l 1241 \7\. That a similar result should hold 
in the non-relativistic setting such as quantum spin systems was long expected and taken for granted 
by theoretical physicists [26j. In [10], Hastings proposed to use Lieb-Robinson bounds to obtain 
such a result and a complete proof was recently given in |2H [T2] . 

As a final application of these locality bounds, we describe a new proof of the Lieb-Schultz- 
Mattis theorem, see [101 122j . in Section [5j These results can be traced back to [10] where Hastings 
introduced a new way to construct and analyze variational states for low-lying excitations of gapped 
Hamiltonians. He developed a notion of a quasi- adiabatic evolution [11] which he then used to 
present a multi-dimensional analogue of the celebrated Lieb-Schultz-Mattis theorem [15]. Such a 
theorem is applicable, for example, to the standard spin-1/2, anti-ferromagnetic Heisenberg model 
and yields an upper bound on the first excited state of order c(log-L)/L for systems of size L. His 
arguments rely on Lieb-Robinson bounds and the Exponential Clustering Theorem in an essential 
way, and we have recently obtained a rigorous proof of this result which holds in a rather general 
setting, see [22] . 

We expect that the ideas currently emerging from recent applications of Lieb-Robinson bounds 
will continue to lead to interesting new results for quantum spin systems in the near future. 

2 Lieb-Robinson Bounds 

The first proof of locality bounds in the context of quantum spin systems appeared in 1972 in a 
paper by Lieb and Robinson [16j. They proved a bound on the group velocity corresponding to 
the dynamics generated by a variety of short range Hamiltonians. In a series of works [21] . [12j . 
and [T9], these estimates have been generalized, and the proof we provide below, see Theorem l2.lt 
illustrates many of the new insights which have recently been developed. 

The result stated in Theorem 12.11 below differs from that which may be found in [16] in 
two important ways. First, the new proof does not require the use of the Fourier transform, and 
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therefore, it extends to models defined on sets without an underlying lattice structure. These results 
may be of interest to those who wish to study quantum spin systems in the context of quasi-crystals 
or in the study of circuits for quantum computation. Second, and most importantly, the constants 
which appear in our bound do not depend on the dimensions of the underlying, single-site Hilbert 
spaces. This opens up the possibility of applying them to models with an infinite-dimensional 
Hilbert space, such as lattice oscillators |20j . 

The basic set up in this theory concerns quantum spins systems, in particular, a finite or infinite 
number of spins labeled hy x €z V. A finite dimensional Hilbert space Tix is assigned to each site 
X £ V. These may represent the spin of an electron, photon, or an atom. In other contexts, 
these states may represent the ground state and first exited state of an atom or a molecule. More 
abstractly, these systems may, for example, model a collection qubits, the basic units of quantum 
information theory and quantum computation. 

If the set V is finite, the Hilbert space of states is given by Tiy = ^xev'^^- -^^^ each spin x, 
the observables are the complex rix x Ux matrices, M„^, where Ux = dim(TCx)- In this context, the 
algebra of observables for the whole system is Ay = 'S)xev -^"^ ■ 

The locality results we wish to describe pertain to observables with finite support. Here, the 
support of an observable is understood as follows. U X C V, we write Ax = 'S>xex ^^n^- 
identifying A S Ax with ^ (g) 11 E Ay, we have that Ax C Ay- The support of an observable 
A £ Av in the minimal set X CV for which A = A' ^ 1 with A' G Ax- 

For infinite V, the algebra of observables is the completion of the algebra of local observables 
given by 

Av= \J -^x 

XCV 

where the union is over all finite X C V. 

To describe the models we wish to investigate, we must first define interactions, local Hamilto- 
nians, and the corresponding dynamics. An interaction is a map <I> from the set of subsets of V to 
Av with the property that ^{X) G Ax and ^{X) = ^{X)* for all finite X cV. A quantum spin 
model is defined by a family of local Hamiltonians, parametrized by finite subsets A gV, given by 

Ht=Yl HX)- (1) 

XcA 

For notational convenience, we will often drop the dependence of on ^. The dynamics, or 
time evolution, generated by a quantum spin model is the one-parameter group of automorphisms, 
{T^}teR, defined by 

T{'{A) = e'^^^Ae-'^^^, A £ Aa, (2) 

which is always well defined for finite sets A. In the context of infinite systems, a boundedness 
condition on the interaction is required in order for the finite-volume dynamics to converge to a 
strongly continuous one-parameter group of automorphisms on Ay- 

The locality results we prove in Theorem 12.11 are valid for a large class of interactions. To 
describe this class precisely, we first put a condition on the set V, which is relevant only in the 
event that V is infinite. We assume that the set V is equipped with a metric d and that there 
exists a non-increasing function F : [0, oo) — > (0, oo) for which: 
i) F is uniformly integrable over V, i.e., 

||F|| = sup < oo, (3) 

and 
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ii) F satisfies 

C= sup V^W^'f/'f'^" <oc. (4) 

x.yGV^ F{d{x,y)) 

Given such a set V, it is easy to see that for any F which satisfies i) and ii) above, then the 
family of functions Fa, for a > 0, defined by 

F„(x)=e-»^F(x), (5) 

also satisfy i) and ii) with ||Fa|| < ||F|| and Ca < C. In this context, we define the set Ba{V) to be 
those interactions ^ on V which satisfy 

^ \mx)\\ 

ma = sup > — — — < OO. (6 

The Lieb-Robinson bounds we will prove are valid for all <I> G Ba{V) with a > 0. Simply stated, 
these results correspond to estimates of the form 

\\[Ti'{A),B]\\ < c(^,S)e-^('^(^'^)~''*l*l), (7) 

where A and B are local observables, t/^(-) is the time evolution corresponding to a finite volume 
Hamiltonian generated by an interaction <I> G Ba{V), and d{A,B) is the distance between the 
supports of A and B. What is crucial in these estimates is that the constants c{A,B), /i, and 

are independent of the volume A C 1^ on which rj^(-) is defined. Physically, the constant 
corresponds to a bound on the velocity at which the dynamics, generated by ^, propagates through 
the system. 

Intuitively, it is clear that the spread of the interactions through the system should depend on 
the surface area of the support of the local observable A, typically denoted by X, not it's volume. 
To make this explicit in our bound, we will denote the surface of a set X, regarded as a subset of 
A C F, by 

Sa{X) = {Z C a : Z n X / and Z n XV 0} • (8) 
Here we will use the notation S{X) = Sv{X), and define 

mux; X) = [ "P^^^ " (9) 

[ otherwise, 

a quantity corresponding to the interaction terms across the surface of X. Comparing the local 
quantity appearing in ([9]) with the norm on the full interaction given by ([6|), one trivially has that 

||$|U(x;X) < mUxdM^), (10) 

where we have used xy to denote the characteristic function of a set y C A, and we have introduced 
the ^-boundary of a set X, written 9<j,X, given by 

a$X = {x G X : 3Z G S{X) with x G Z and $(Z) / } . (11) 

The Lieb-Robinson bound may be stated as follows. 
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Theorem 2.1 (Lieb-Robinson Bound) Let a > and take ^ G Ba{V). For any finite set 
A C denote by t/^(-) the time evolution corresponding to the local Hamiltonian 

H^=Y^ HZ). (12) 

ZcA 

Given any pair of local observables A G Ax o-nd B G Ay with X,Y <Z K and d{X, Y) > 0, one may 
estimate 

\mA\B]\\ < [e^mWA _ D^i^xXh (13) 

for any t G K. Here the function Da{X,Y) is given by 



DaiX,Y)=mm 



(14) 



A number of comments are useful in interpreting this theorem. First, we note that if X and Y 
have a non-empty intersection, then the argument provided below produces an analogous bound 
with the factor e'^ H*"" I*' — 1 replaced by H*"" l*L In the case of empty intersection and 
for small values of \t\, (jl3p is a better and sometimes useful estimates than the obvious bound 
||[rt(A),S]|| < 2P||||5||, valid for all t G M. 

Next, if $ G Ba{V) for some a > 0, then one has the trivial estimate that 

DaiX,Y) < \\Fo\\\ma min(|a$X|,|5$y|) e-'"^^^'^K (15) 

Clearly then, we have that 

llr A^zi^ Hill ^ \\B\\ IIFqII . ^a\d(X,Y)-'-^^\t\\ 

lirt (^),^]|| < ^ mm {\d<s>X\, \d<s>Y\) e L « J, (16) 

which corresponds to a bound on the velocity of propagation given by 

< inf 3^MSL, (17) 

a>0 a 



Next, we observe that for fixed local observables A and B, the bounds above, p3j) and (|T6]) . 
are independent of the volume A C given that A contains the supports of both A and B. 
Furthermore, we note that these bounds only require that one of the observables has finite support; 
in particular, if \X\ < oo and d{X, Y) > 0, then the bounds are valid irrespective of the support of 
B. 

Proof of Theorem \2.1[ To prove (fT3ll , we will provide an estimate on the quantity 

AeAz ll^ll 

where i? is a fixed observable with support in Y, and the subset Z C A we regard as arbitrary. 
Introduce the function 

f{t)=[Ti^{r'',{A)),B], (19) 

where A and B are as in the statement of the theorem. Due to the strict locality of the Hamiltonian 
Hx, as defined e.g. in (jl2p and the fact that the observable A G Ax, we have that supp(r^((^)) C X 
for all i G M. It is straight forward to verify that 

f'it)=i E [rHHZ)),f{t)]-i E [rtHr^tiA)),[THHZ)),B]]. (20) 
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(21) 



Since the first term in (j20p is norm preserving, we find that 

\\[THr^M)),B]\\ <\\[A,B]\\+2\\A\\ ^ f \\[rHHZ)), B]\\ ds. 

The inequahty (I2ip and the fact that ||r^((^)|| = II^H together imply that 

CBiX;t) < CBiX;0) + 2 J] mZ)\\ CB{Z;s)ds. (22) 

zas{x) 

It is clear from the definition, see (jlSp . that for any finite Z C A, 

Cb{Z;{)) < 2\\B\\6y{Z) (23) 

where Sy{Z) = OifZny = and 6y{Z) = 1 otherwise. Using this fact, iteration of ()22p yields 
that ^ 

Cb(X,0 < 2||S|| Vffl^an, (24) 

where for n > 1, 

n 

«n= 5^ 5^ ••• E ^y(^n)n 11^(^^)11- (25) 

For an interaction $ S ;Ba(y), one may estimate that 

«i ^ EE E 11^(^)11 ^ EEii^ii'^(^;^)^'^(^(^'y))- (26) 

x&X y&Y zes(x): x&X y£Y 

x,yeZ 

In addition, 

«2 < EEE E 11^(^1)11 E ii'^(^2)ii 

xeX yeY zeA ZieS(X): Z2&S{Z^): 
x,z£Zi z,yGZ2 

x&X yGY zeA ZieS(x): 

x,z£Zi 

^ W'^w- E E E (^(^' ^)) 2/)) 11*^11'^ ^) 

xex yeY zeA 

< \maCaY,Y.m\a{x;X)Fa{d{x,y)), (27) 

xeX y£Y 

where we have used dH) for the final inequality. With analogous arguments, one finds that for all 
n > 1, 

an < miaCaT'' J] ^ || ||a (x; X) F„ (d(x, y)) . (28) 

a-eX y& 

Inserting ([25]) into (|24p we see that 

Cb(X,0 < ^{{^ (e^W^^^I^I-l) j;||<I>|U(x;X)F,(d(x,y)), (29) 

xGX y£Y 

from which (I13p immediately follows. □ 
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3 Quasi-locality of the dynamics 

The Lieb-Robinson bounds of Theorem 12.11 imply that the dynamics of quantum lattice systems, 
those generated by short range interactions, are quasi-local in the sense that the diameter of the 
support of any evolved local observable does not grow faster than linearly with time, up to an 
arbitrarily small error. There are at least two interesting ways to give precise meaning to this quasi- 
locality property of the dynamics. In the first, one shows that the time-evolved observable can be 
well-approximated in norm by one with a strictly local support. This is achieved by the quantum 
version of integrating out the variables in the complement of a ball with a radius proportional to 
time. In the second, we show that to compute the dynamics up to a time t > 0, one can replace 
the Hamiltonian with a local Hamiltonian supported in a ball of radius proportional to t. Clearly, 
the net result is the same: the support of observables evolved with approximate dynamics remains 
contained in the ball where the local Hamiltonian is supported. After presenting the details of these 
two approaches, we conclude this section with a few interesting applications that immediately follow 
from quasi-locality. 

As in the previous section, we will work with an interaction <I> G Ba{V) with a > 0. For the 
purpose of the discussion below, we will consider a finite subset A C ^ and restrict our attention to 
the dynamics r^^(-) generated by the finite volume Hamiltonian H\ as defined in ([T]); our bounds 
will be independent of the volume A. For any X C A we will denote by = A\X. 

In the first approach to obtaining a local approximation with support contained in A" C A, one 
takes the normalized partial trace over the Hilbert space associated with X^. In order to estimate 
the norm difference it is convenient to calculate the partial trace as an integral over the group of 
unitaries [5]. Given an arbitrary observable A £ Aa and a set X C A, define 



where IA{X'^) denotes the group of unitary operators over the Hilbert space Tix'^ and /i is the 
associated, normalized Haar measure. It is easy to see that for any A G Aa, the quantity {A)x': 
has been localized to X in the sense that {A)x<^ £ Ax- Moreover, the difference may be written in 
terms of a commutator, i.e. as 



To localize the dynamics, let A S Ax and fix e > 0. Based on our estimates in Theorem 12. H 
we approximate the support of t/^(v1) with a time-dependent ball 




(30) 





(31) 




(32) 



For any unitary U G U{Bt{e, A)^), we clearly have that 




(33) 



a 



and therefore, using ([3T]) above and our bound (fT6|) . we immediately conclude that 




2\\A\\ \d^X\ 

C'a 




< 



^o||e- 



(34) 
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The tolerance e > can be chosen to optimize estimates. 

In the second approach, one shows that only those terms in the Hamiltonian supported in 
BT{e,A) contribute significantly to the time evolution of A up to time T [19j. Again, we consider 
the finite volume dynamics applied to a local observable A E Ax- Fix e > 0, T > 0, and consider 
the ball BT{e,A) as defined in (j32p above. The estimate 



(A) 



(35) 



valid for all |t| < T, readily follows from the results in p9|. In fact, the proof of ()35p uses the 
following basic estimate, see e.g. Lemma 3.3 in [19] , 



Lemma 3.1 Let $ G Ba{V) with a > and take a finite subset A C V. If the finite volume 
Hamiltonian is written as the sum of two self-adjoint operators, //a = +H^^\ and for i = 1,2, 



(i) 

denotes the dynamics corresponding to H)^ , then the following estimate is valid: 



rtiA) - Ti'\A) 



< 



ds, 



(36) 



for any observable A and t G 



To apply Lemma [3. II in the context discussed above, we write the local Hamiltonian as the sum 
of two terms: 



Ha 



(37) 



ZfSAtBTie,A}) 



ZCA: 

ZeSAiBT{e,A)) 



Recall that for any X C A, we defined S\{X), the surface across X in A, with equation dS]). Drop- 
ping the surface terms comprising H^^ above, decouples the dynamics, i.e., t^^\A) = tI^'^^^''^\A), 
and we find that for any \t\ < T, 



tHA) - T, 



BT{e,A) 



(A) 



< 



E 



ZeSAtBrieA)) 



Brie, A) 



(A) 



ds. 



(38) 



For each of the terms on the right hand side above, the Lieb-Robinson estimates imply that 



Bri^A) 



(A) 



< 



2||$(Z)|| ||yl|| 

C'a 



(39) 



as 



Careful consideration of the combinatorics of these sums, exactly as done in [TU], yields 
claimed. 

It is interesting to ask about locality bounds for complex times. For one-dimensional systems 
with finite range interactions, Araki proved that the support does not grow faster than an expo- 
nential in \z\, where z is the complex time [31 Theorem 4.2]. In other words, the complex time 
evolution 



-R,R] 



(A) 



AzHa 



Ae 



can be approximated by ri ^'^\A) with small error, uniformly in i?, as long as log \z\ < cN, for a 
suitable constant c. For stochastic dynamics of classical particle systems, good locality bounds that 
are very similar to Theorem 12.11 are known [171 Theorem 4.20]. Since these classical models with 
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stochastic dynamics can be thought of as quantum systems, generated by particular Hamiltonians 
and evolving with purely imaginary times, this may indicate that the general result, which allows 
for no more than logarithmic growth of should not be considered the final word on this issue 
in regards to specific models of interest. 

The quasi-locality formulation of the Lieb-Robinson bounds makes it easy to derive bounds on 
double and higher order commutators of the form 

[TtM)ArtAB),Tt,{C]]. 

Such commutators remain small in norm as long as ti,t2i and are such that the intersection of 
the time-dependent "quasi-supports" remains empty. That is whenever ti,i2)^3 £ 1^ are such that 

Bt,{e,A)r\Bt,{e,B)r\Bt,{e,C) = (k. 

Another immediate application of the quasi-locality property is a bound on the rate at which 
spatial correlations can be created by the dynamics starting from a product state. In p[9], we 
proved the following. 

Theorem 3.2 Let a > 0, $ G BaiV), and take Vt to he a normalized product state. Given X, y C A 
with d{X,Y) > and local observables A £ Ax and B £ Ay, the estimate 

\{Tt{AB)) - {n{A)) {Tt{B))\ < A\\A\\ \\B\\ \\F\\ {\d^X\ + \d^Y\) e""'^^^'^) (40) 

is valid for all t G M. Here, for any observable A, the expectation value of A in the state Vt is 
denoted by (A) = {i},Ai}), and the function 

Gait) = 9^±}M ||ci>||, e^ll^lU^^N ds. (41) 

Jo 

4 Exponential Clustering 

As a second application of these locality bounds, we will present a proof of the Exponential Clus- 
tering Theorem, see Theorem 14.11 below, which improves on the estimates found in [2T|, see also 
|12| . The proof of exponential clustering demonstrates that models with a spectral gap above the 
ground state energy necessarily exhibit exponential decay of spatial correlations in their ground 
state. Such results have recently appeared in a variety of contexts, e.g. [5T1 [T^]. Using valence 
bond states, as is done in [18], one can easily construct gapless models with exponentially decaying 
ground state correlations indicating that, in general, there is no converse to Theorem 14.11 

In the finite volume, the notion of a gapped Hamiltonian is clear. If the system is infinite, we 
express the gap condition in terms of the limiting dynamics, the existence of which is guaranteed 
by the Lieb-Robinson bounds as discussed above, by considering a representation of the system on 
a Hilbert space TC. This means that there is a representation vr : Ay — > S{Ti.), and a self-adjoint 
operator H on TC such that 

TTiniA)) = e''"TT{A)e~''^, A e Av (42) 

We will assume that the representative operator H is non-negative and that there exists a vector 
€ Tl for which HQ = 0. We say that the system has a spectral gap in the representation if there 
exists 5 > such that spec{H) n (0, 6) = 0. In this case, the spectral gap, 7, is defined by 

7 = sup{5 > I spec{H) n (0, 6) = 0}. (43) 

Let Pq denote the orthogonal projection onto kerH. From now on, we will work in this represen- 
tation and simply write A instead of 7r{A). 
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Theorem 4.1 (Exponential Clustering) Let a > and take <I> G Ba{V)- Suppose that the 
dynamics corresponding to ^ onV can he represented by a Hamiltonian H with a gap 7 > above 
the ground state energy, as described above. Let Q be a normalized ground state vector for H ; 
i.e. satisfy HVt = with ||0|| = 1. Then, there exists a constant ^ > such that for any local 
observables A £ Ax and B £ Ay with X,Y CV and d{X, Y) > satisfying PqBQ = PoB*Q = 0, 
the bound 



\{n,Anh{B)n)\ < C{A,B,^)e 
is valid for all non-negative b satisfying < 67 < 2fj,d{X,Y). One may take 

07 



mWaCa + 7' 



as well as a constant 



C{A,B,-f) 



\A\\ \\B\ 



1 



+ 



2||Fn 



fid{X,Y) TlCa 



mill (|5$X|, \d^Y\ 



(44) 



(45) 



(46) 



Note that in the case of a non-degenerate ground state, the condition on B is equivalent to 
{Q, BQ) = 0. In this case, the theorem with 6 = becomes 



I {n, ABn) - {n, aq) {n, Bn) \ < c{a, b, 7) e 



-^id(X,Y) 



which is the standard (equal-time) correlation function. For small 6 > 0, the estimate 
viewed as a perturbation of (j47p . Moreover, for 6 > large, there is a trivial bound 



Ath,{B)Q)\ < \\A\\ \\B\\ e 



-by 



(47) 
can be 

(48) 



Proof of Theorem \4.1\ 

We will follow very closely the proof which appears in [21j and refer to it whenever convenient. 
Consider the function / given by 



f{z) = {n,AT,{B)n) 



f 



e'''' d{A*n,PEBn) 



(49) 



where we have used the spectral theorem and the fact that B projects off the ground state. It is 
clear that the function / defined in (I49p is analytic in the upper half plane and has a continuous 
(and bounded) boundary value on the real axis. The quantity we wish to bound corresponds to 
f{ib) for 6 > 0. The case 6 = will follow by a limiting argument. 

Since the boundary value of / on M is continuous, one may show by a limiting argument that 
for any T > 6, 



fiib) 



1 

27ri 



Ft 



ib 



-dz. 



(50) 



where Tt is the semi-circular contour from —T to T (on the real axis) into the upper half plane. 
As is shown in [2T], the fact that the Hamiltonian is gapped, see also ()49p . implies that the integral 
over the circular part of the contour vanishes in the limit as T — > oo, and therefore we have the 
bound 



\{VL,Anb{B)Q)\ 



\f{ib)\ < limsup 

T^oo 



1 

27ri 



T 



-T 



fit) 
t - ib 



dt 



(51) 
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The proof is then complete once we estimate this integral over the real line. While this inequality 
is true for any value of 6 > 0, to get the desired estimate, we will have to choose 6 > sufficiently 
small; see the comments following (j64p below. 

Let a > 0. We will choose this free parameter later. Observe that one may write 



fit) 



-at-' 



(52) 



Clearly then, the integral we wish to bound, the one on the right hand side of ([5T|) above, can be 
estimated by 



1 

27ri 



T 



fit)e 



T 



t - ib 



-dt 



+ e" 



1 
27d 



fit) 



T 



t - ib 



dt 



(53) 



We will bound the absolute value of each of the integrals appearing in (j53p separately; the prefactor 



will be an additional damping made explicit by the choice of a. 



To bound the first integral appearing in (j53p . we further divide the integrand into two terms. 
Note that 

/(t)e-"*' = (J],rt(i?)^J7)e-°*' + (f], [A rt(i?)] 0)6""*' . (54) 
By the spectral theorem, we have that 



1 

27ri 



{n,TtiB)AQ)e 
t - ib 



-dt 



CO Y 
2TTi 



T ^—itE^—at^ 



t - ib 



dtd{PEB*n,An), 



(55) 



where we have used now that i?*f2 is also orthogonal to the ground state. Applying Lemma |4.2 
stated below, to the inner integral above, we have that 



1 



lim 



t - ib 



dt 



1 



, , e e 4" dw < -e 4a 

2^/^./o ~ 2 



(56) 



where for the final inequality above we used that E > j > 0, a > 0, and 6 > 0. From this we easily 
conclude that 



lim sup 



1 

27ri 



T 



{n,TtiB)An)e- 

t - ib 



-dt 



\\A\\ \\B\\ _7i 
< e 4" . 



For the integral corresponding to the second term in (I54p . it is easy to see that 



(57) 



27ri 



{n,[A,TtiB)] n) e~ "*^ 
t - ib 



-dt 



< 



2lT 



WiAniB) 



-at' 



dt, 



(58) 



where we have taken advantage of the fact that 6 > 0. To complete our estimate, we will introduce 
another free parameter s > 0. Here we use the Lieb-Robinson bound, Theorem l2.1l for times |t| < s 
and a basic norm estimate otherwise. The result is that the right hand side of (|58p is bounded from 
above by 

^"^""^"-D,(X,y)fe2ll*ll'^^'^^-l) + Ml^e— ^ (59) 



This completes the bound of the first integral appearing in (I53p . 
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Using again the spectral theorem, the second integral in (j53p may be written as 



T e 



itE I „ab'^ 



-T 



t - ib 



dtd{A*n,PEBn) 



(60) 



As is described in detail in [21], we find that for E >j and a chosen such that 7 > 2ab, 



1 fT e^iS 
lim / 

T^oo 27ri t — ib 



dt < —e 4" , 



(61) 



which produces an estimate (one analogous to the bound in (f57|) above) for (f60]) . 

All of our estimates above combine to demonstrate that the right hand side of (j5ip is bounded 

by 



„ .„ „„„ r 2Da{X,Y) / 2ii$ii r • 

\\A\\\\B\\ e~i^ + — 111 II V. 11^11°^°^ 

if a satisfies 7 > 2ab. The choice a = 7/2s yields: 

7£ 

B e-^ 



1 + 



1 



1 I 2 ^ 2D,(x,y) J.ii^ii^c^. 



(62) 



(63) 



As is demonstrated in (|15p . Da(X, y) decays exponentially as e '*'^('''^'^). In this case, if we choose 
s to be the solution of the equation 



s{2\maCa + 7/2) = adiX,Y), 



(64) 



then we have proven the result. Notice that we have chosen a in terms of s, which is defined 
independently of b, thus the condition 7 > 2ab will be satisfied for sufficiently small 6 > 0. □ 
In our proof of the Exponential Clustering Theorem above, we used several times the following 
useful fact, a proof of which appears in |21j . 



(65) 



Lemma 4.2 Let E e R, a > 0, and z G C+ = {z G C : Im[z] > 0}. One has that 



1 e'^t^-at^ 

lim / dt 

T-»oo 2TTi J_x t — Z 



1 



2v7ra 



e e 



Moreover, the convergence is uniform for z £ C 



5 The Lieb-Schultz-Mattis Theorem 

As a final application of these locality bounds, in particular both Theorem 12.11 and Theorem 14.11 
we were recently able to provide a rigorous proof the Lieb-Schultz-Mattis theorem, see |23) . which 
is valid in arbitrary dimensions. In this section, we will discuss this result and outline the ideas 
which motivate our proof. 
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5.1 The Result and Some Words on the Proof 



The classical Lieb-Schultz-Mattis Theorem (LSM), [15], concerns the spin-1/2, anti- ferromagnetic 
Heisenberg chain. This model is defined through a family of Hamiltonians H]^^ acting on the Hilbert 
space = (2)x'G[i L] '^ifh the form 



L-l 



Hl = J2Sx-S,+i. (66) 

x=l 

Here, for each integer x E [1,-^^], the spin vectors Sx have components 

5^ = •••®11®5^ ®11®---® 1, i = l,2,3, (67) 

where is the corresponding spin-1/2 (Pauli) matrix acting on the x-th factor of Wji^/,]. The LSM 
Theorem may be stated as follows. 

Theorem 5.1 (LSM, 1961) If the ground state of Hl is unique, then the gap in energy between 
the ground state and the first excited state is bounded by C/L. 



A further result by Lieb and Mattis in 1966, [14], verified that under certain conditions, for 
example when L is even, the main assumption in the LSM Theorem, specifically the uniqueness of 
the ground state, is indeed satisfied. Almost twenty years later, the LSM Theorem was generalized 
to encompass a variety of other one (and quasi-one) dimensional models by Affleck and Lieb in pQ . 
In particular, this result applies to those chains of even length with spins having arbitrary half- 
integer magnitude. Note: Here and in the rest of this section the term half-integer, or half-integral, 
refers to one-half of a positive odd integer, i.e., an element of the set N -|- 1/2. 

For models to which the LSM Theorem applies, one expects that the excitation spectrum cor- 
responding to the thermodynamic limit has no gap above the ground state energy. It is interesting 
to note that the predictions of Haldane ^ suggest that such a result is rather sensitive to the type 
of interaction terms. In fact, the spin-1, anti-ferromagnetic Heisenberg chain is predicted to have 
a robust gap above the ground state energy in the thermodynamic limit. 

In a work of 2004, see [10], Hastings argued that a higher dimensional analogue of the LSM 
Theorem could be proven using the improved locality bounds which have recently been established. 
We will now summarize these ideas and indicate how they may be implemented to demonstrate a 
rigorous proof of this theorem. 

The multi-dimensional LSM Theorem, stated as Theorem 1 5 . 2 1 b elow . is valid for a large class of 
models; a detailed proof of this is contained in [23]. For simplicity of presentation in this review 
article, we will restrict our attention to the spin-1/2, Heisenberg anti-ferromagnet, however, our 
general assumptions are discussed in Subsection I5.2.4[ In u dimensions, the model of interest is 
defined on subsets V;, C Z'^ in analogy to (j66p above, i.e. one considers Hamiltonians 



Hl= ^ Sx ■ Sy 



x,yeVi^: 

\x-y\=l 

acting on the Hilbert space = ^xgVl state the new result. 

Theorem 5.2 // the ground state of Hi is non- degenerate, then the gap, jl, above the ground 
state energy satisfies 

7. < C'-^. (69) 
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The logarthmic correction which appears in Theorem 15.21 in contrast to the original result of 
Theorem 15. 1^ seems to be an inevitable consequence of the locality bounds we incorporate in our 
proof. It is an interesting open question to determine whether or not there is a class of models, in 
dimensions > 1, for which one can prove such a bound without the logarithmic correction. 

In essense, Theorem 15.21 is proven using a variational argument. Letting ipQ denote the unique, 
normalized ground state, we know that for any normalized vector ipi with | (-(/'o j V'l ) I 1; the bound 



is always valid. Here we have dropped the dependence of all quantities on the length scale L. Prom 
this perspective, there are only three steps necessary to prove the desired result. First, we must 
construct a normalized trial state ipi, as indicated above. Next, we must estimate the difference in 
the energy corresponding to ipi and Eq. Lastly, we must ensure that the inner product |(^0)V'i)l 
remains sufficiently small. 

This method of proof is complicated by the fact that the ground state is virtually unknown, and 
therefore the means by which one should construct a trial state is not apriori clear. Inspiration for 
the construction of our variational state comes from the work of Hastings, again see [10], in which 
he proposes to consider the ground state of a modified Hamiltonian, Hg, where the interactions in 
a given hyperplane have been twisted by an angle of 6; more on this below. The ground state of 
this modified Hamiltonian may be regarded as the solution of a specific differential equation, in 
the variable 9, whose initial condition corresponds to the unique ground state whose existence we 
assumed. The solution of Hastings' differential equation is ameanable to analysis, in particular, one 
can apply both the Lieb-Robinson bounds and the Exponential Clustering Theorem to provide the 
desired energy and orthogonality estimates mentioned above. One may recall that the clustering 
bounds, as in Theorem 14.11 provide estimates which themselves depend on the size of gap 'Jl- For 
this reason, the argument proceeds by way of contradiction. In fact, by assuming that there exists 
a sufficiently large constant C for which the gap satisfies -fL > C\og{L)/L for large enough L, we 
construct a trial state whose energy eventually violates this bound. 

5.2 A More Detailed Outline of the Proof 
5.2.1 Constructing the Trial State 

In our proof, we use the fact that the Hamiltonians we consider are assumed to have at least one 
direction of translation invariance. We incorporate this into our notation by considering finite 
subsets Vl C Ij^ of the form V/, = [1, L] x where we have isolated a particular direction, which 
we will often refer to as the horizontal direction, and perpendicular sets C with cardinality 
|V^| < CL^~^. For the orthogonality result, we will also need to assume that |V^| is odd, see 
Subsection 15.2.31 The trial state is constructed from a perturbation of the Hamiltonian Hi defined 
by "twisting" certain interaction terms. A twist in the hyperplane situated at a site m £ [1,L] is 
defined by replacing all interaction terms Sx ■ Sy in (j68p corresponding to horizontal bonds with 
X = [m, v), y = {m + 1, v), and some v E Vf^ by terms of the form 




(70) 




(71) 



for some G R. A doubly twisted Heisenberg Hamiltonian is then given by 




(72) 



\x-y\=l 



14 




Figure 1: Plot of the three lowest eigenvalues of Hq q, the singly twisted Heisenberg Hamiltonian 
for a ring of 8 spins. The magnitude of spins is S = 1/2 in the plot on the left and S" = 1 on the 
right. 



where 

6, if X = (m, v),y = {m + l,v) for some v G V^, 
ea:y=le', if X = im + L/2,v),y = {m + l + L/2,v) ioi V £Vf^, (73) 
0, otherwise. 

Here we have taken periodic boundary conditions in the horizontal direction. 

It is interesting to note the behavior of the lowest eigenvalues of the singly twisted Heisenberg 
Hamiltonian Hq^q for a simple spin ring with an even number of spins. The behavior depends in 
an interesting way on the magnitude of the spins. When the spins are half-integer, the two lowest 
eigenvalues cross at = n. In contrast, when the spins are integer, they remain non-degenerate. 
The quasi-adiabatic evolution is a device designed to construct a continuous path from the ground 
state of Hq q to the first excited state of i?27r,o which, of course, are both identical to the unperturbed 
Hamiltonian. 

One of Hastings' crucial insights in [TO] is that, for the half-integer spin case, the first excited 
state can be obtained by applying a "quasi-adiabatic evolution" to the ground state; here 9 is 
the evolution parameter. More concretely, let il)o{9,9') and Eq{9,9') denote the ground state and 
ground state energy of the doubly twisted Hamiltonian Hq^q', respectively. It is easy to see that along 
the path 9' = —9, the Hamiltonian Hg _g is unitarily equivalent to the unperturbed Hamiltonian 
Hl = Hofi. By differentiating the eigenvalue equation Hq_qiIjq{9,—9) = Eo{9,—9)'ipo{9,—9) and 
using that dgEQ{9, —9) = 0, one obtains 

deMG, -0) = - „ ^ „ [deHe,-e\MO, -0)- (74) 
ne-e — ^0 

Formally, equation (j74p may be re-written using the Heisenberg dynamics as follows, 

deM0,-0) = B{9)M0,-0), (75) 
where the operator B{9) is defined by 

POD 

B{9) = - Tu{deHe,-e)Po{9,-9)dt. (76) 
Jo 
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Here we have denoted by Tt{-) the dynamics generated by the Hamiltonian Hg^^o, and Po{9, —9) is 
the corresponding spectral projection onto the ground state. Due to the gap assumption, (I76p is 
well-defined. 

Equation (175p captures the evolution of the ground state of a doubly twisted Hamiltonian along 
the path 9' = —9 where the effect of the first twist is canceled by the second. The trial state is 
obtained by using an approximation of the differential equation in (I75p to describe the ground state 
of the singly twisted Hamiltonian Hg Q. Introduce the parameters a > and T > 0. To approximate 
the imaginary time evolution corresponding to an arbitrary Hamiltonian H, with dynamics Ts{-), 
at time t > 0, we define 

for any observable A. In general, an anti-Hermitian operator of the form 

Bo^^riA, H) = - fiA^iit, H) - A^{it, H)*] dt, (78) 
Jo 

will be used in place of ([76]) . 

Note that the observable which is evolved in equation (j76p . dgHg^-g, contains terms of two 
types. The first are localized around the twist of angle 9 which correspond to horizontal bonds 
{x,y) of the type x = {m,v), y = {m + l,f), and some v G ^l- '^^^ second are similar, yet 
localized around the twist of angle —9. We group these two types of terms together and write 
dgHg^^Q = diHg _g — d2Hg _g to simplify notation. By linearity, the operator 

Ba,T {dgHg^^g, Hg^^g) = Ba,T {dlHg^^g, Hg^^g) - B^^t {92^^8-9, Ho-o) ■ (79) 

Given a sufficiently large gap 7/,, there is a choice of the parameters a and T for which the ground 
state ipo{9, —9) is well approximated by the solution of the differential equation (|75p with B{9) 
replaced by Ba^T {dgHg _g,Hg _g). Hastings' clever proposal is that one may also approximate the 
ground state of the singly twisted Hamiltonian, V'o(^) 0), by evolving with just the first term on the 
right hand side of ()79p above. More concretely, consider the operator Ba^T{9) defined by setting 
Ba,T{9) = Ba^xidiHg^-g, Hg^-g) and solve the Hastings' Equation given by 

dg^Pa,T{(^) = B^,T{9)lPa,T{G) (80) 

with initial condition V'o,t(0) = ^o{0,0), i.e., the ground state of the unperturbed Hamiltonian. 
Under the choice of parameters a = Jl/L and T = L/2, the resulting variational state V'l = 
■i/'a,r(27r) may be estimated in such a way that Theorem 15.21 follows. We note that a particularly 
nice feature of the differential equation (|80p is that the operator Bo,^t{9) is anti-Hermitian, and 
therefore the solution remains normalized for all 9 G [0, 27r]. Hastings interprets the solution of (|80p 
as a quasi- adiahatic evolution of the ground state V'o = ''Ao(0, 0). 

5.2.2 Locality and the Trial State 

One key technical lemma, which we use repeatedly in all of the estimates that follow, concerns the 
local evolution of the solution to ([80]) . Consider a sub- volume A£,(m) C Vl localized around the 
twist of angle 9, for example, take Ai(m) to be of the form [m — {L/A — 2),m + {L/A — 2)] x 
which is slightly less than half of the system. Let Pa,T{9) and pq{9, —9) denote the density matrices 
corresponding to the states 'ipa,T{9) and ipo{9, —9). 
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Lemma 5.3 Suppose there exists a constant c > such that L^l > c and choose a = 'Jl/L and 
T = L/2. Then, there exists constants C > and k > such that 

sup \\Try^\AUm) [Pa,m " Poi9, -0)]\\, < CL^^e'^'^^K 

e6[0,27r] 

Since the proof proceeds by way of contradiction, the assumption L^fi > c is part of the 
argument. Moreover, it is easy to produce a bound on 7^ from above that is independent of 
the length scale L. Lemma 15.31 demonstrates that if L^l is sufficiently large, then the effect of 
ignoring the second twist in the definition of ipa,T{G) is negligible when one restricts their attention 
to observables localized in AL{m). The proof of this lemma uses both the Lieb-Robinson bound. 
Theorem 12. H and the Exponential Clustering result. Theorem 14. II 

5.2.3 The Estimates 

The Energy Estimate: To estimate the energy of the trial state, we consider the function 

E{e) = (Va,T(^), He,-ei'a,Tm. (81) 

Due to the initial condition used to define ipa,Tid), we know that E{0) = Eq the ground state 
energy, and since H2tt-2tt = Hi, E{2-k) corresponds to the energy of the trial state ipi = ipa,T{'^'^)- 
The main idea here is to use the locality property of the trial state, i.e. Lemma [5.3| and the unitary 
equivalence of the Hamiltonians Hg _g to obtain an estimate on the derivative of this function which 
is uniform for 9 £ [0, 2tt]. Explicitly, we can prove the following result. 

Theorem 5.4 Suppose there exists a constant c > such that L^l > c and choose a = Jl/L and 
T = L/2. Then, there exists constants C > and k > such that 

KVi,i^LVi)-i5^o| <CL3-ie-'=^^^ 

The Orthogonality Estimate: As we mentioned before, Hastings' quasi-adiabatic evolution is 
norm preserving. In particular, we are guaranteed that ||'0i|| = HV'oll = 1- Our argument that 
■01 is sufficiently orthogonal to Vo makes essential use of the fact that the total spin in each 
perpendicular set, V^, is half-integer. In the case of the spin-1/2, anti-ferromagnetic Heisenberg 
model, this corresponds to the assumption that |V^| is odd. We have the following theorem. 

Theorem 5.5 Suppose there exists a constant c > such that L^l ^ c and choose a = 'Jl/L and 
T = L/2. Then, there exists constants C > and k > such that 

KVi.V'o)! <CL2-e-'=^^^ 

To prove this result we observe that, although the ground state ipo{9,—6) of the perturbed 
Hamiltonian is not translation invariant, it is invariant with respect to "twisted" translations. 
In fact, let T be a unitary implementing the translation symmetry in the horizontal direction, 
specifically the direction in which we have imposed periodic boundary conditions, chosen such that 
TipQ = ipQ. This is possible since '00 is the unique ground state and the Hamiltonian is translation 
invariant with respect to T, i.e., T*HlT = Hl. 

Define twisted translations by setting 

Te,e' = TUrrA0)U^+L/2{0') (82) 



17 



where the column rotation, Un{0), apphes the rotation e*^*^^ to all sites for the form x = (n, v), for 
some V € V^. The unitary equivalence of the doubly twisted Hamiltonian to Hi = Hq^q can 

also be expressed in terms of these column rotations 

He,-e = W{erHo,oW{e), (83) 

where 



Wie)= 09 Uni9). (84) 

m<n<m+L/2 

With these definitions, it is easy to see that W{6)*TW{6) = Tq^^q commutes with Hg^^g, and 
therefore, 

Te,-eMO,-0) = MO,-0), (85) 

as we claimed. 

The main idea in the proof of Theorem is to again use the locality properties of the solution 
of Hastings' Equation (j80p to show that 

Te,Otl^a,T{^) ~ V'a,T(^). (86) 

Since the total spin in each V^^ is half-integer, again a consequence of assuming |V^| is odd, the 
column rotation [/m(2vr) = —1. Clearly then, T2n,o = —T, and therefore, Tipi ~ —tpi- As we have 
chosen T so that, TipQ = -00, this implies that ipi is nearly orthogonal to ipQ. 



5.2.4 Generalizations 

As we previously indicated, the proof of the multi-dimensional Lieb-Schultz-Mattis Theorem, 
demonstrated in ^23j, applies to a large variety of models. In what follows below, we will out- 
line a list of assumptions which define a wide class of Hamiltonians for which the LSM Theorem 
remains valid. 

The Basic Set- Up: It is not important for our argument that the underlying sets have a lattice 
structure, in particular the sets Vl need not be subsets of Z"^ . Rather, we need only assume that 
there is, at least, one direction of increase, which previously we labeled the horizontal direction. 
We make this notion concrete by assuming that there exists an increasing sequence of finite sets 
{Vl}l>i, exhausting some infinite set V, which are of the form Vl = [^,L] x V^. We also assume 
a bound on the cardinality of the perpendicular sets of the form \ V^\ < cL°' for some a > 0, and 
it is natural, but not necessary, to take a = — 1. 

The interactions can also be of a general form. We assume that the set V is equipped with 
a metric d and a function F as described in Section [2j To start with, we work with interactions 
<1* G Ba(y) for some a > so that the infinite volume dynamics is well defined from the beginning. 
In contrast to (168p . the more general finite volume Hamiltonians are of the form 

= ^ ^{X) + boundary terms, (87) 
XcVl 

where we will assume periodic boundary conditions in the horizontal direction and arbitrary bound- 
ary conditions in the other directions. 

Assumption I. Our first assumption is that the interaction <I> G Ba{V) is translation invariant 
in the horizontal direction. This is clearly the case for the Heisenberg anti-ferromagnet, and in 
general, it means that for any X C Vl, 

<^>iX + ei) = aimX)), (88) 
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where X + ei is the translation of all points in X by one unit in the horizontal direction and ai(-) is 
the translation automorphism which maps Aj^^ y-±) into Aj^^j^i y±) for all n G Z. Here the column 

sets (n, V^) are defined by (n, V^) = {x ^ Vl : x = {n, v) for some v G V^}- Due to the assumed 
periodicity in the horizontal direction, this translation invariance can be implemented by a unitary 
T e ^Vi, i-e. ^{X + ei) = T*^{X)T for ah X C Vl- This unitary T will depend on the length 
scale L, but we will suppress this in our notation. 

Assumption II. We further assume that the interaction has a finite range i? > in the 
horizontal direction. This assumption is not strictly necessary. It is clear from the estimates in |23] 
that the result remains true even if the interactions are of long range with sufficiently fast decay. 

Assumption III. We assume the interaction $ has rotation invariance about one axis. Again, 
this is clearly the case for the Heisenberg model. In the more abstract setting, we specifically 
assume that for each x ^ V there is a local hermitian matrix, which we will denote by S"^, with 
eigenvalues that are either all integer or all half-integer. These matrices are also required to be 
translates of one another, i.e., for any x £ V, ai{S^) = Rotation invariance for a general 

interaction <I> means that for any X C Vl, 

U*{e)^{X)U{e) = $(X) for all G M, (89) 

where the rotation U{9) is defined by 

U{e) = (5?) e*^^'. (90) 



Assumption IV. We assume that the matrices introduced above, i.e. S^, are uniformly 
bounded in the sense that there exists a positive real number S for which sup^g^ ll'S'^H < S. 
In addition, we must assume an odd parity condition on the spins. Define the parity px of a site 
X G y to be if the eigenvalues of S'^ are all integers and 1/2 if they are all half-integers. The odd 
parity assumption is that 

P(n,.) eN+1/2, (91) 

for all n G Z. For the spin-1/2 Heisenberg model, we satisfied this assumption by taking the 
cardinality of the perpendicular sets, |V^|, to be odd. In general, the sum of the spins over the 
perpendicular set needs to be half- integer to ensure that the column rotations Un{9), as defined 
after equation ([82]) . satisfy [/„(27r) = — 1. As we have seen, such an identity plays a crucial role in 
our argument for orthogonality. 

Assumption V. The ground state of Hi is assumed to be non-degenerate. In this case, it is 
also an eigenvector of the translation T and the rotations U{9) introduced above. We assume that 
the ground state has eigenvalue one for both T and U{9). 

Assumption VI. We assume that there are orthonormal bases of the Hilbert spaces with 
respect to which and ^{X) are real for all x £ V and all finite X C V. This assumption is 
satisfied by the Heisenberg model, and it is an important symmetry which allows us, in general, to 
prove that the ground state eigenvalue is invariant with respect to the doubly twisted Hamiltonian 

The following theorem was proven in [23]. 



Theorem 5.6 Let <I? G Ba{V) for some a > 0. If ^ satisfies assumptions I-VI above, then the gap, 
'jL, above the ground state energy of Hi satisfies 

7. < CM^. (92) 
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